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Abstract 



J> \ We find and solve a large class of integrable dynamical systems 

t" which includes Calogero- Sutherland models and various novel gener- 

alizations thereof. In general they describe ./V interacting particles 
■<^j- \ moving on a circle and coupled to an arbitrary number, m, of su(iV) 

spin degrees of freedom with interactions which depend on arbitrary 
■ real parameters Xj, j = 1,2, ... ,m. We derive these models from 

SU(iV) Yang-Mills gauge theory coupled to non-dynamic matter and 
on spacetime which is a cylinder. This relation to gauge theories is 
used to prove integrability, to construct conservation laws, and solve 
these models. 

o 

Integrable models have always played a central role in classical and quan- 
tum mechanics. Most prominent examples, like the Kepler problem, are sys- 
tems with few (< 3) degrees of freedom. An important exception is a class of 
integrable A-particle models associated with the names Calogero, Moser and 
Sutherland 0] (for review see ||). These are models for identical particles 
moving on one dimensional space and interacting via certain repulsive two- 
body potentials v (r). A well-known example is v(r) oc g 2 / sin 2 (gr) (which 
includes v(r) oc 1/r 2 in the limit g — > 0), and we refer to the correspond- 
ing model as CS model. It is known that these models allow for interesting 
generalizations which also have dynamic spin degrees of freedom 0, ||. The 
CS model and its generalizations have recently received quite some attention 
in different contexts. Here we only mention their relation to gauged matrix 
models |J and gauge theories on a cylinder [0 which will be relevant for us.0 

a the latter relation is implicit already in earlier work; see e.g. |q] 
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In this article we find and solve a new class of integrable systems con- 
taining the CS model and their spin-generalizations as limiting cases. 
Our method is to extend and exploit the relation of the CS model to gauge 
theories on a cylinder, as will be explained in detail below. For simplicity 
our discussion here is restricted to classical models, and we only consider a 
special type of gauge theory. Our method is simple, and it should be pos- 
sible to generalize it in several different directions: We conjecture that the 
corresponding quantum models are also integrable, and that our method to 
prove integrability should apply to the quantum case, too. (In this context it 
is worth noting that the quantum-analog of the gauge theory we consider is 
closely related to QCD on a cylinder and in the limit where the masses of the 
quarks becomes infinite; a good starting point to the literature on this is Ref. 
H). Moreover, it would be interesting to extend our method to other gauge 
theories (e.g. with supersymmetry and/or more general types of matter fields 
than the one considered here) and thus try to find and solve other integrable 
models. 

The models we solve are given by a Hamiltonian 

N / ct\ 2 i JV m i N m 

h = E %r + ~ E E - /)pf pT + n E E <wOT 

0=1 a,p=\ j,k=l Z a=l j,k=l 

where q a and p a are particle coordinates and momenta with the usual Poisson 
brackets {q a ,p 13 } = S al3 etc., and p°f = pj a are complex valued su(iV)-spins, 
i.e. Ea=oPr = 0and 

{pf, Pf'} = igtoS* (S**pf - 5^ P f) . (2) 

(The other Poisson brackets vanish.) The interaction potentials are given by 

where 

x jk = (xj - x k ) 27T (4) 
with S2tt '■= s — 27m for the integer n such that — tt < s — 27m < tt, and 



The parameters g (real positive), N and m (positive integers) are arbitrary, 
and — 7r<£i<£2<...< x m < 7r. Furthermore we have the following 
constraint on the possible initial conditions, 

in 

EpT = Q V «- ( 6 ) 

Note that Vjk{r) = Vjk(-r) = r), which implies that the Hamiltonian 
is real. Moreover, since 

v jk (r + ^) = e- m ^- x ^v jk (r), (7) 

the Hamiltonian is invariant under the following transformations, 

q a ^q a + — , p a ^p a , pf - p f e -^ a -^) (8) 
9 

for all integers n a . Thus these models describe particles moving on a circle 
of length 1/g and interacting with a potential whose strength depends on 
dynamic spins. We note that the particles repel each other, and we have a 
further natural restriction on phase space, 

q a ^qP Va^f3. (9) 

The main result of this article is a proof of integrability and the explicit 
solution of all these models. 

Next the relation of our particle-spin models to gauge theories is dis- 
cussed. It is known that the CS model can be obtained from a gauged one 
dimensional matrix model 0. More recently a relation to gauge theories in 
1+1 dimensions was pointed out In this article we explore this relation 
further and use it to find and solve new integrable models. We present a sim- 
ple argument that SU(iV) Yang-Mills gauge theory on a cylinder coupled to 
certain non-dynamic matter is equivalent to a model of interacting particles 
and spins. We then show that this equivalence can be used as a powerful tool 
to analyze and solve these models: integrability is manifest, the construction 
of conservation laws trivial, and a simple solution method is obtained by ex- 
ploiting gauge invariance, i.e. the possibility to change from the Weyl gauge^ 

b Ao and A\ will be defined further below 
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Aq = to what we call the diagonal Coulomb gauge, i.e. the condition that 
the spatial component of the Yang-Mills field, Ai(t,x), is independent of x 
and diagonal in color space, 

A 1 (t,x) = Q(t)=diag(q\t),q 2 (t),...,q N (t)) . (10) 

This is due to the fact that the gauge theory model in the Weyl gauge is 
free and can be solved trivially, whereas in the diagonal Coulomb gauge the 
time evolution equations are non-linear and, in a special case, equal to the 
Hamilton equations of the model given by Eq. ([[]). To be more specific: We 
restrict ourselves to gauge theory models with matter fields localized at a 
finite number of points Xj, j = 1, 2, ... m for simplicity (see Eq. ( p5|) below). 
We find that the dynamics of such a model in the diagonal Coulomb gauge 
is governed by the equations of motion X = {X, H}, X = q a ,p a , Pj , which 
follow from the Hamiltonian Eq. ([!]) and the Poisson brackets {■,■} given 
above. These observations allows us to derive the full solution of the initial 
value problem for all these models generalizing the known solution of the CS 
model 0. 

We now consider 1+1 dimensional Yang-Mills theory, i.e. the differential 
equations X)«=o,i[-^mj F^ v ] = J u where = + igA^, A^ are the Yang- 
Mills fields, J M the matter currents, g the Yang-Mills coupling strength, = 
d/dx^ 1 , and p, v = 0, 1. Spacetime is a cylinder i.e. t = x° G R is time, and 
x = x 1 G [— 7r, 7r] (= circle). Moreover,^ F^ u = D u ]/ig, and our metric 
tensor is diag(l, —1). As gauge group we take SU(iV) in the fundamental 
representation.^ 

We restrict ourselves to non-dynamic matter, i.e. J 1 = 0, and J° = p. We 
denote p as charge. Note that we have to impose [D , p] = for consistency. 
Setting E := F i, we can write these equations as follows, 

d Q A 1 = E + d 1 A Q + ig[A x ,A Q \ (11) 

d E + tg[A ,E] = (12) 

d p + ig[A , p] = (13) 

d 1 E + ig[A 1 ,E]=p. (14) 

Eq. (^) is called Gauss' law and is a constraint on possible initial data for 
c [a, b] := ab — ba 

d i.e. Afj,, E and p are functions with values in the traceless, complex N x N matrices 
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the system of time evolution equations (|TT|)-(|r3|). We now exploit gauge in- 
variance: Eqs. (|ll])-(0) are obviously invariant under gauge transformations 

A„ -> U^A^U + -U~ l dJJ 

w 

E -> U~ X EU 

p -> C/^V^ (15) 

where {7 = U(t,x) is an arbitrary differentiable SU(iV)-valued function on 
spacetime. To eliminate the gauge degrees of freedom one has to fix a gauge. 
One convenient choice is the Weyl gauge A (t,x) = 0.| Then the Eqs. (pTT|)— 
can be solved trivially: 



E(t,x) = E(0,x), 
A x (t,x) = A 1 (0,x) + E(p,x)t, 
p(t,x) = p(0,x) (16) 

with the initial data E(0,x), Ai(0,x) and p(0,x) satisfying the Gauss' law 



Eq. Q14|) (note that our solution Eq. ([H]) satisfies Eq. (|14]) for all t if it satisfies 
it for t = 0). 

As mentioned, E, A^ and p are functions with values in the traceless NxN 
matrices. In the following we write the matrix elements of M = E, A^ or p as 
M a/3 , a, (3 = 1, 2, . . . N. Note that, since J2a=o M ° a is zero > the independent 
components are M a/3 for a^, and M aa - M a+1 > a+1 for a = 1, 2, ... iV - 1. 

We now show that one can also impose the diagonal Coulomb gauge (]10D , 
i.e. for each (generic) Yang-Mills configurations Ai(t, x) one can find a gauge 
transformation U such that A\ = U~ 1 AiU + U~ 1 d\U /ig is a diagonal matrix 
Q independent of x [JTOj . For that we construct such a U explicitly. We first 
note that a solution to the equation d\S + igA\S = with S(t, —n) = 1 is 
the parallel transporter 

S(t,x) =Pexp (-^ £ dyAifayj) (17) 

where P exp is the path ordered exponential. Note that S(t, x) is not a gauge 
transformation since it is not periodic in x (its values at x — —tt and tt are 

e i.e. to consider the model in terms of the gauge transformed fields on the r.h.s. of 
Eq. (|l5|), which by abuse of notation we denote by the same symbol, and with a gauge 
transformation U which is a solution of doll + igAoU = 0. 
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different in general). To construct a gauge transformation, we introduce the 
SU(iV)-matrix V(t) diagonalizing the SU(iV)-matrix S(t,ir), 

vity^it, 7r)v(t) = e - i92nQ{t) (is) 

for some diagonal matrix Q(t). This implies that 

U{t, x) = S{t, x)V(t)e l9(x+n)Q(t) (19) 

is periodic in x, and it satisfies d\U + igA\U = igUQ equivalent to A± = Q. 
Moreover, if Ax(t,x) is a generic differentiable map on spacetime, then Q(t) 
and V(t) can be chosen to be differentiable in t [O, and U(t,x) Eq. (|Tj]) 
is indeed a differentiable function on space-time i.e. a gauge transformation. 
'Generic' here means that the latter is only true if q a (t) 7^ q^(t) for all t 
and a//J since otherwise discontinuous functions V(t) can occur [|TTJ. Due 
to Eq. @, gauge field configurations Ax(t,x) where this condition fails are 
irrelevant for us. Note that our discussion here implies that the q a (t) can 
be obtained as eigenvalues of the Wilson line S(t,n). This observation will 
allow us to determine the explicit solution of the Hamilton eqs. following 
from Eqs. ([!]) and (|3]). 

We now determine the time evolution equations for the variables q a de- 
fined in Eq. (0). We use Fourier transformation, 

E al3 (t,n)= r dxe- inx E a P(t,x), neZ (20) 

J — IT 

and similarly for A and p. Then Eq. ( |TTD gives 

d q a (t)=p a (t) = Ea y> 0) . (21) 

Z7T 

Note that this and the following equations all are consistent with J2a=i = 
J2a=i P a = (this corresponds to translation invariance of the mechanical 
system defined in Eq. ([!])). The time evolution of the p a follows from Eq. 

(0), 



in N , 

doP a (t) = --^ E E (Af(t,n)E^(t,-n)- 

\ Z7T ) neZ 13=1 
f3^a 



-E a P(t,n)A% a (t,-n)) . (22) 
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The r.h.s. of this equation can be evaluated using Eqs. ( ITT] ) and ( fH| ) 

- i (n + s[g Q (t) - g^)]) Af (t, n) = E a ?(t, n) 

i (n + s[g Q (t) - <A*)]) £ a/3 (t, n) = p a/3 (t, n) (23) 
(note that this holds true even for a = j3 if n 7^ 0). Inserting this we get 

Note that up to now no specific choice for the charges was made. To 
proceed, we restrict ourselves to charges of the following form for simplicity, 

m 

p^(t,*) = X>f (25) 
which describe matter localized at the points Xj, as mentioned above. Then 

m 

p^( tj n) = J2pf(t)e- inx \ (26) 

k=l 



and we can then write Eq. (|33J) as 

Q N m 

9opa V = "TwTT E E pfW/^CtKC^W - /W) (27) 

eg l r J 0=1 j,fe=i 

with 

e m(xj-a: fe ) 

%W = E7^r— 1^ (28) 
V (2vr) 2 (n + #r) 2 

where the summation is over all n G Z. This equals Eq. (|3|), as can be seen 
by a simple computation using the identity 

pins / 2 \ 

E 7~T^ = ( —27^ + ^ " * W ) • (29) 

Note also that 



neZ ( n + r ) \sin (7rr) 



e ins sL , , 7T 2 



Et7 = T-W + T' (3°) 
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Eqs. (^1|) and (|27| ) are precisely the Hamilton equations q a = {q a ,H} and 
p a = {p a -i H}. We are left to determine the time evolution of the p°^ . From 
Eq. (|T3p we get 

JV 

do P f(t) = -ig £ (K\t, xj)pf(t) - pT(t)Af(t, Xj j) . (31) 

7=1 

Moreover, we can compute A^it.Xj) = ^J2 n Ao l3 (t,n)e mXj from Eqs. ([23]) 
and Eq. (|26|) . Note that Eq. ( p3|) also determines AQ a (t,n) if n is non-zero, 
and we can set A^it, 0) = 0.[] Then a simple computation gives 

m 

Af(t, x s ) = 2vr £ « iJfc (g°(f) - <f(t))p?{t) (32) 
fc=i 



with Ujfc(r) given by Eq. (j|) and Ujfc(0) = by Eq. (0). (We used Eqs. (|28|), 
( |29|) and (|30l)). Note that for a = [3 the summation in Eq. (|^) is restricted 
to the non-zero integers n, which implies Vjk(0) = Cjk-) With that Eq. (PT|) 
becomes equal to the Hamilton eq. p a ^ = {p^,H} following from Eqs. ([!]) 
and (0). 

We finally note that the n = components of Gauss' law Eq. ( |14"D for 
a = (3 reads 

p aa (t,0) = Va. (33) 

This is a consistency requirement. Our arguments above show that this 
condition is fulfilled for p in Eq. (|25| ) if and only if Eq. @ holds for all t, 



which is true if it holds for t = 0. 

We now solve the Eqs. (|l|), (g7|) and (^TJ) (g2j) with the initial conditions 

<f(0) = C P a (0)=Po, pf(0)=/$. (34) 

Our discussion above implies that we can obtain the solution q a (t) of this 
initial value problem by solving the Yang-Mills Eqs. ([ITD-(p^) with the initial 
conditions 

r7T dx 



„ fir fix 

Af(t = 0,x) = 6 aP <fi, / — E aa (t = 0,x 

J-TT 271 



Po 

N 



f ^(t = 0,x)='£p$5(x-x J ). (35) 



•'A non-zero Aq a (t) = is irrelevant since it can be removed by a gauge transformation 
compatible with the DCG. 
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We first have to determine E a ^(t = 0,x) for a ^ (3 from Gauss' law Eq. 
(|i4|) . The solution A±(t, x) of the gauge theory in the Weyl gauge A = 
is then given in Eq. flip]). To obtain the q a (t), we only need to evaluate 
the corresponding parallel transporter S(t, 7r) Eq. (|17|): as discussed, the 
eigenvalues of S(t,n) are equal to e~ 2mgqa ^\ Moreover, 

Pj (t) = U&x^pwU&Xj) (36) 

with U(t,x) given by Eq. (jl9|). Here and in the following we use an obvious 
matrix notation. 

For t = we can write Gauss' law Eq. flT4| ) as follows^] 

d 1 (e i9Qox E{0, x)e~ i9Qox ) = e i9Qox p{0, x)e~ i9Qox (37) 

with p(0, x) = jyjLi Pj,o$( x ~ ^i) an d 

Qo = diag(ql,q%,...qo) . (38) 

Since p(0,x) = except for x = Xj, we obtain E(0,x) = e~ l9< ^° x Bje l9< ^° x , 
where -B, is some constant matrix, for Xj < x < Xj + i, j = 0, ...,m, 
x = —Ti and x m+ i = n. To determine the matrices Bj we integrate Eq. 
( |37|) from Xj — + to Xj + + . This gives the recursion relations Bj — 
Bj_i = e % ^ x i pjQe~ ig Q° X:i ', and the condition E(0,—tt) = E(0,tt) implies 
e 2igQ n j^^ e -2igQ TT _ Q m Putting this together and using the second relation 
in Eq. ( |35l) we obtain after a straightforward calculation, 

i=\ i=\ Z7r e=i j 
(l-ME— a ( 39 ) 

£=1 2«sin(p7r[gg ! - q'o}) 



with sgn(x) = 1 for x > and —1 for a; < 0. With that, Eq. ( |16|) gives 
ar) = e - 49Qo:E (Q + Bjt) e igQox for < x < x j+1 . This is the solution of 
the Yang-Mills equations in the Weyl gauge. 

We now solve d\S(t,x) +igAi(t,x)S(t,x) = which is equivalent to 

d\S(t, x) + igBjtS(t, x) = for Xj < x < Xj + \ (40) 



9 Note that the same argument applies for times t > 
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for S(t,x) = e i9QoX S{t,x). This implies S(t,x) = e- i9B ^ x - x ^S(t,Xj) for 
Xj < x < Xj + i, thus 

S(t,x j+1 ) = < "'^'•' ( '■■'S{l..r J ) = ...= 

X • • • e -i9 B ot{xi+n) e ~igQon /^-q 

where we used S(t, —n) = 1. Especially (for j = m), 

^ _ e -igQoir e -igB m t(ir-x m ) _ _ _ 
xe -«A-ii(i».-im-i) . . . g-igBotix! +?r) g-igQoT (42) 

We thus obtain our main 

Result: The Hamiltonian equations of the the dynamical system defined in 
Eqs. ar e given in Eqs. fl2l[ ) and ft3l\). The solutions of these 



equations with the initial conditions Eq. P^ j can be obtained from the eigen- 
values of the matrix S(t, 7r) given in Eqs. and ftSQ ) according to Eq. 
(fj^). Moreover, pj(t) is given by Eq. ( |5^[ j with U(t,Xj) defined in Eq. ( |7^ j 
and S(t,Xj) in Eq. 



Note that for m = 1, the dynamics of the spin and the particles decouple, 
and our result reduces to the known solution of the CS model; see e.g. 0. 

It is now also easy to construct conservation laws for our dynamical sys- 
tems: Eq. (|T6|) implies that tr[E(t, x) n ], where tr is the N x N matrix trace, 
is time independent for all — n < x < n and all positive integers n. Since 
these quantities are gauge independent, they are time independent also in 
the diagonal Coulomb gauge. In this latter gauge, we can evaluate E(t,x) 
as above and obtain E(t,x) = e~ tg ®^ x ' Bjit^e 19 ®^ for Xj < x < Xj + \ where 
Q(t) and Bj(t) are as in Eqs. (|38|) and ([$9|) but with q% , and p°g replaced 
by q a (t), p a {t), and (t), i.e. the solution of the initial value problem which 
we solved above. Using cyclicity of the trace, we conclude that tr[Bj(t) n ] for 
an arbitrary positive integer n and j = 1, . . . m are time independent: Each 
of them is a conservation law. For m = 1 these are the known conservation 
laws for the CS model [[]]. It is also worth noting the corresponding Lax-type 
equations d Bj(t) + ig[Mj(t), Bj(t)} = where 

Mj(t) = e igQ{t)x > A (t, Xj )t "' Qu)r ' - XjP{t) (43) 
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which are obtained from eq. fll^) setting x = Xj and using doQ(t) = P(t) : = 
diag(pi(t),p 2 {t), ...,p N (t)) (A {t,Xj) is given by Eq. (|3|)). 

As discussed above, the models we find generalize the CS models with 
the interaction potential v (r) oc g 2 / sin 2 (gr) which describe particle moving 
on a circle of length 1/g. There is a an integrable CS-type model of particles 
moving on the real line and interacting with a potential v (r) oc g 2 / sinh 2 (gr), 
see e.g. ||. The sinh- model and its solution can be obtained from the sin- 
model and its solution by replacing q a — > iq a and p a — > ip a [[3]. This re- 
placement in our Hamiltonian, Eq. (H) (together with H — > —H), leads to 
spin-generalizations of the sinh-model. It is natural to conjecture that this 
very replacement allows to obtain the solution of the latter from our solution 
of the former model. 
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